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ABSTRACT 

We introduce several useful utilities in development 
for the creation and analysis of real wildland fire 
simulations using WRF and SFIRE. These utili- 
ties exist as standalone programs and scripts as 
well as extensions to other well known software. 
Python web scrapers automate the process of down- 
loading and preprocessing atmospheric and surface 
data from common sources. Other scripts simplify 
the domain setup by creating parameter files au- 
tomatically. Integration with Google Earth allows 
users to explore the simulation in a 3D environment 
along with real surface imagery. Postprocessing 
scripts provide the user with a number of output 
data formats compatible with many commonly used 
visualization suites allowing for the creation of 
high quality 3D renderings. As a whole, these im- 
provements build toward a unified web application 
that brings a sophisticated wildland fire modeling 
environment to scientists and users alike. 

Index Terms — Data preprocessing, Data anal- 
ysis, Geophysics computing, Client-server systems 

1. INTRODUCTION 

The weather research and forecasting model (WRF) 
is a full-featured mesoscale weather model used 
as a research tool and operational forecasting (TJ. 
SFIRE is an add-on wildland fire forecasting model 
coupled with WRF J2j . An earlier version is avail- 
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able in WRF release as WRF-FireQ Because WRF 
is widely used in the research community, there are 
a large number of standard tools available for creat- 
ing and visualizing its simulation outputs. However, 
SFIRE acts on domains that are significantly higher 
resolution than the typical WRF simulation. In ad- 
dition, SFIRE uses surface variables on a refined 
subgrid that is different from the standard atmo- 
spheric surface variables. These subgrid variables 
and the higher resolution modeling typically re- 
quires special processing for input and visualization 
that is not supported by most commonly used tools. 
For this reason, we have begun to develop new 
software that mimic and extend standard utilities 
for data processing and visualization specifically 
tailored for use with SFIRE. 

2. MODEL INITIALIZATION AND DATA 
PREPROCESSING 

As computers become more powerful, there is in- 
creasing interest in simulating smaller scale phe- 
nomena using WRF than is typically associated 
with mesoscale weather forecasting. In particular, 
a typical SFIRE simulation occurs at mesh resolu- 
tions on the order of 10 m or less. Even the highest 
resolution surface data provided with WPS is sev- 
eral hundred times coarser. For these fine scale 
domains, the ability to import custom datasets into 
WPS is essential for the initialization of a realis- 
tic simulation Resources such as the USGS's 
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seamless data server that provide open access to 
high quality surface data are generating a large in- 
terest in initializing WRF simulations using custom 
datasets. WPS provides a mechanism for these 
datasets through a simple binary file format that is 
described in the WRF technical documentation [4]; 
however, there is no standard API or GIS software 
capable of writing to it. Users who lack sufficient 
technical knowledge are currently unable to process 
this data into WRF. 

Prior work related to WRF-Fire has lead to 
small utilities that are able to convert standard GIS 
GeoTIFF files into Geogrid's binary file format. 
The GeoTIFF specification is commonly used and 
can be provided as output from a wide range of GIS 
applications. TopoGrabbe^jis a Python application 
based on this work that is capable of download- 
ing and converting topological data automatically. 
More recently, modifications to Geogrid have been 
written allowing it to read GeoTIFF files directly, 
and extensions to WPS standard programs have 
made it possible to ingest GeoTIFF data directly 
from the USGS into a WRF simulation (5]. The 
implemented changes to Geogrid allows a user to 
import a USGS GeoTIFF dataset directly into the 
WRF workflow without the (often) difficult and 
error prone procedure of converting the data. The 
GeoTIFF interface is compiled in as an optional 
component and is useful for both experienced GIS 
users with complicated workflows and those who 
just want to replace a single data source. The im- 
plementation allows overriding erroneous metadata 
without modification to the source images, and effi- 
cient access to the data. 

With tools in place to automatically retrieve and 
integrate known data sources into a simulation, it 
is now possible to automate nearly all of the steps 
necessary to initialize a simulation for WRF and 
SFIRE. A command-line client has been written us- 
ing Python to generate all required parameter and 
input files for a simulation given user defined do- 
main and fire ignition specifications. This script is 
currently being integrated with a web application 
based on Google Maps, providing users with a fully 
graphical environment for designing a simulation on 
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real surface imagery. This client application is cou- 
pled with a server located on a high performance 
compute cluster that runs the simulation on demand 
and returns results back to the client as they are pro- 
cessed. 

3. MODEL EVALUATION AND DATA 
POSTPROCESSING 

By default, WRF generates NetCDF binary files 
containing the raw simulation output. These files 
contain a large number of variables that are useful 
in analyzing the simulation output; however, they 
are often too large to be transferred to the user's 
computer for visualization. Several utilities have 
been written that allow users to visualize the sim- 
ulation without downloading the full output. For a 
quick overview of the fire propagation, it is possible 
to generate a KMZ file containing a raster image 
of any surface variable or a georeferenced polygon 
outlining the burning region at a given timestep. 
This KMZ file can be opened in a number of GIS 
applications including Google Earth to be overlaid 
on a rendering of the earth's surface. Google Earth 
allows the user to play back the images in the KMZ 
files as an animation with a timestamp indicating 
the simulation time. The web application in devel- 
opment uses this capability to display the results of 
the running simulation on top of Google Maps. 

There are a number of options for rendering a 
3D visualization of the simulation. VAPOR is an 
application developed at NCAR designed for visu- 
alization of standard WRF atmospheric simulations. 
The postprocessing utility wrf2vdf takes a stan- 
dard WRF NetCDF output file and generates a more 
compact vdf file containing only the fields the user 
wishes to investigate. VAPOR features a very fast 
and easy to use GUI for generation of high qual- 
ity images; however, it does not support the refined 
grids used by the fire variables. When using VA- 
POR with SFIRE, users are limited to low resolution 
versions of fire variables, which can result in blocky 
images. 

Other popular 3D visualization applications in- 
clude Paraview, MayaVI, and VisTrails. These ap- 
plications all use the Visualization Toolkit (VTK) 



as a backend and do not have any built in capabil- 
ity for reading WRF output files. SFIRE contains 
a Python script wrf 2vtk that can generate VTK 
compatible files from WRF outputs at full resolu- 
tion. This script is much like VAPOR's wrf 2vdf 
in that it converts the data to a more compact repre- 
sentation and compresses the files for transfer to a 
local computer. 

4. TOWARD A FIRE FORECASTING WEB 
PORTAL 

Automated scripting of the model setup procedure 
has lead to the development of a web portal provid- 
ing users with a simple interface for running SFIRE 
simulations using real data. The user registers with 
the site with a username and password. Upon log- 
ging in, he or she is presented with a Google Maps 
interface with instructions to click on a point on the 
map to create a new ignition. When clicked, a dia- 
log box appears allowing the user to customize the 
parameters of the simulation. The user is notified 
by email on completion of the simulation when the 
user can return to the web portal to view the output 
of this and all prior simulations. Currently, the out- 
put of the simulation displays as an animation the 
heat flux into the atmosphere from the fire. In the 
future, we plan to add the ability to visualize other 
surface variables such as fuel types, winds, and fire 
danger ratings. 

The web frontend is built on top of a Google 
Maps interface using its javascript interface. The 
core web application including user and content 
management is served by Django. The backend 
component is executed at a remote site on a com- 
pute cluster accessible via ssh from the web server. 
The web application passes a string in the form of 
simple keyword=value pairs specifying parameters 
of the fire simulation. These parameters include 
the ignition location and time chosen by the user, 
domain size and resolution, and simulation run time 
among others. The web server executes an appli- 
cation script on the compute server providing the 
simulation parameters and a session ID. The com- 
pute server initializes and executes the simulation 
while a post processing daemon watches for new 



output files. Whenever a new output file is detected, 
a visualization script is executed generating a KMZ 
file, which is then uploaded to the web server to be 
displayed to the user. 

5. CONCLUSION AND FUTURE WORK 

The utilities presented here are designed to enhance 
the usability of SFIRE for the average user by elimi- 
nating the need for complex data conversion or an in 
depth understanding of the underlying code. They 
are leading to a unified application where a user can 
design a simulation through a series of graphical di- 
alogs, spawn the simulation on a remote compute 
cluster, and visualize the output of the running sim- 
ulation locally. The development of such an appli- 
cation would provide an invaluable tool for research 
and education. 
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Fig. 1. SFIRE web portal results page. Along with a standard Google Maps interface, the results portal 
contains animation controls so the user can view the development of the fire front. Prior simulations are 
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John Michalakes, Syed Rizvi, Xin Zhang, 
Jonathan D. Beezley, Janice L. Coen, and 
Jan Mandel, "ARW version 3 modeling 
system user's guide," Mesoscale & Mis- 
croscale Meteorology Division, National 
Center for Atmospheric Research, July 
2010, |http : / / www . mmm . ucar . eduTI 
|wrf /users /docs /user_guide_V3/ 
|ARWUsersGuideV3 . pdf | 

[5] Jonathan D. Beezley, A. K. Kochanski, V. Y. 
Kondratenko, and J. Mandel, "Integrating 
high-resolution static data into WRF for real 
fire simulations," Paper 6.3, Ninth Sym- 
posium on Fire and Forest Meteorology, 
Palm Springs, October 2011, 2011, |http:| 



1/ / ams . conf ex . com/ ams/ 9FIRE/ 
|webprogram/Paper!92276 . html[ 
retrieved December 201 1. 



